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In the first section, I will try to convey a sense of the variety of observational 
inputs that tell us about the existence and the spatial distribution of dark 
matter in the universe. In the second section, I will briefly review the four 
main dark matter candidates, taking note of each candidate's status in the 
world of particle physics, its production in the early universe, its effect upon 
large scale structure formation and the means by which it may be detected. 
Section 3 concerns the energy spectrum of (cold) dark matter particles on 
earth as may be observed some day in a direct detection experiment. It 
is a brief account of work done in collaboration with J. Ipser and, more 
recently, with I. Tkachev and Y. Wang. 
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Ref. [1] is a list of works which I have consulted in preparing this rather cur- 
sory overview, and which the reader should turn to for more complete and in-depth 
information. 

I. DARK MATTER OBSERVATIONS 

In 1932, Oort^ studied the motion of galactic disk stars in the vertical direction, 
i.e., perpendicular to the disk. By applying a version of the virial theorem to the 
distribution of vertical star velocities, he obtained an estimate of the density of the 
galactic disk in the solar neighborhood: 

Pdi^^l^-lO-^gr/cm 3 . (1.1) 

On the other hand, if one adds up the densities of all the matter "seen" in stars and 
interstellar gas, plus what is expected from stellar remnants, mainly white dwarfs, one 
finds considerably less than the dynamical estimate of Eq. (1.1), of order half thereof. 
So there is dark matter in the galactic disk. Because this dark matter is in the disk 
rather than the halo, we expect it to be dissipative, which means in all probability that 
it is baryonic dark matter. 

In 1933, Zwicky^l used measurements of the line-of-sight velocities of galaxies 
in the Coma cluster to estimate the mass of that cluster using the virial theorem. The 
result he obtained in this way is approximately 400 times the mass inferred by counting 
the number of galaxies in the cluster and assigning to each a mass (~ lO n M ) typical 
for the luminous part of a spiral galaxy. The masses of the luminous parts of some 
spiral galaxies had already been determined by measuring their rotation curves up to 

distances from their centers of order their disk radii. SmithW obtained a similar result 
for the Virgo cluster. 

In 1973, Ostriker and Peebles t 5 l pointed out that the tendency of galactic disks 
to be unstable towards a large-scale bar mode can be cured by assuming the existence of 
a spherical halo of dark matter with mass within the disk radius of order the disk mass 
(6 • 1O 1O M0 for our galaxy). Thus the galactic mass within a sphere with radius equal 
to the disk radius would be roughly half in the disk and half in an unseen spherical 
halo. Of course, by Birkhoff 's theorem, the argument does not say anything about halo 
matter outside the disk radius. 

During the seventies, the rotation curves of spiral galaxies were measured^ 
over much larger distances than before, in many cases extending the rotation curve to 
distances several times the disk radius. In all cases, the rotation velocity was found to 
be constant (i.e., independent of radius r) or slightly rising, up to the last measured 
point. Balancing centrifugal and gravitational forces, one has 

GM{r) = v rot (rf 

where M(r) is the mass interior to r and v ro t(r) is the rotation velocity at r. If only 
luminous matter were contributing to the galactic mass, we would have M(r) ~ constant 

and hence v ro t(r) ~ T2, for r > disk radius. Instead, the data show v ro t(r) ~ constant 
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there and hence M(r) ~ r. The implication is that there is a halo of dark matter whose 
density Pdm( r ) ~ a ^ large r. The halo distribution is usually modeled by the function 

PdmW = " 7772 ( L3 ) 

^ \a) 

where a is called the core radius. For our own galaxy, v ro t — 220km/s, a ~ few kpc, 
and : 

Pdm(r Q ) * \ ■ 10- 24 gr/cm 3 (1.4) 

where r© ~ 8.5 kpc is our distance to the galactic center. The estimate (1.4) of the local 

dark halo density is based upon models^ of the galactic mass distribution developed 
in the early 80's by Bahcall and Soneira, and Caldwell and Ostriker. However, the 
discovery of an abundance of microlensing events in the direction of the galactic bulge 
has stimulated a lot of recent work on the galactic mass distribution and this will likely 
result in a more precise determination of the galactic halo parameters. 

In 1972, J. Einasto, A. Kaasik and E. Saart 8 ] studied 105 pairs of galaxies, 
the members of each pair being close on the sky and assumed to be gravitationally 
bound to each other. They compared the distance between each pair to its relative 
velocity to obtain an estimate of its reduced inner mass. Of course, only the line-of- 
sight velocities and the angular projections of distance onto the sky are measured and 
therefore an average over many pairs must be performed to try and eliminate the effects 
of projection and ignorance of orbit eccentricities. At any rate, these authors find that 
the galactic mass increases with distance, approximately linearly, up to masses of order 

1O 13 M . This study and others^ imply that galactic halos extend very far out. I do 
not know of anything that contradicts the assumption that galactic halos extend all the 
way to radii of order 1-2 Mpc, where the halo density, falling off as becomes equal 
to the average intergalactic dark matter density. 

There are a number of methods to estimate the average dark matter density on 
scales larger than the typical intergalactic or intercluster distance (> 10 Mpc). Density 
perturbations on such large scales are still in the linear regime of their growth by gravi- 
tational instability. Let us describe a particular method. If a region has an overdensity 
Sp in excess of the average density p, neighboring galaxies will have an excess gravi- 
tational attraction towards that region and consequently deviate from perfect Hubble 
flow. One writes: 

v = Hor + v p (1.5) 

where Hq is the Hubble expansion rate, r is the position relative to the center of an 

overdensity and v p is called the peculiar velocity. It is found ^ that in the linear regime 
around a single, spherically symmetric overdensity 

v p = -# fV 6 ^ , (1.6) 
P 

where 

0=^ = ^. (1.7) 

Pcrit 3Hq 
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p cr it is the critical density for closing the universe. The ft dependence on the RHS of 
Eq. (1.6) is a close fit to the actual ft dependence for vanishing cosmological constant. 
Eq. (1.6) affords a way to determine ft by measuring peculiar velocities v p and com- 
paring them with observed overdensities However, ^ cannot be measured directly. 
What can be done is count galaxies, measuring their average density uq and local over- 
densities Suq. Unfortunately, the relationship between ^ and is not known. It is 
parametrized by a fudge factor 6, called the "bias parameter" : 

^ = 6^. (1.8) 
n G p 

So the method of peculiar velocities actually measures ^j— . A number of authors t 10 ! 
have analyzed galaxy distributions in this way with the result: 

ft - 6 

— — = 1±0.3. (1.9) 
b 

Most attempts to determine the bias parameter from first principles yield b > 1. 
Eq. (1.9) suggests then that, when measured on the largest scales, the value of the 
density parameter ft is consistent with a critically closed universe (ft = 1). On the 
other hand, the measurements on these large scales are very imprecise. 

Of course, ft = 1 is strongly favored on 'theoretical' grounds. AO^l universe 
will deviate from = 1 more and more as time goes on, in pretty much the same way as 
a pencil standing nearly vertically on its point will fall over. For a universe to stay near 
ft = 1 for a long time, it has to be extraordinarily close to ft = 1 to start with. This 
problem of initial conditions for our universe is called the flatness (or age) problem. It 
may be neatly solved by assuming that there is, at very early times, a brief epoch during 
which the energy density is dominated by vacuum energy density and, as a result, the 
universe expands at an exponential rate. After this "inflation" , ft = 1 with tremendous 

precision. The inflationary cosmology has many other attractive features as wellJ 11 ! So 
there are compelling reasons to believe that ft = 1. Whether observations support this 
prejudice is not obvious, although it seems fair to say that they are in rough agreement 
with it. Luminous matter contributes ft| UTO — 3.1CP^ to 6.10 - ^. Dark matter in galactic 
halos and in clusters of galaxies contributes Vt ga i ~ 0.02 to 0.2. Finally, as Eq. (1.9) 
and the results of other observations on the largest scales studied suggest, there may be 
enough dark matter not associated with galaxies or clusters of galaxies (in voids, say) 
to yield ft = 1. 

The success of nucleosynthesis t 12 l in producing the primordial abundances of 
light elements requires that the contribution ft^ of baryons satisfies: 

o.on < o.on/r 2 <n B < o.oi9/i~ 2 < 0.12 (1.10) 

where h parametrizes the present Hubble expansion rate 

ff = 100* • . (1.11) 

s Mpc 
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Measurements of Hq are in the range of 0.4 < h < 1. Since f^ um < 0.006, Eq. (1.10) 
implies that some baryons are dark. Recall that there is dark matter associated with 
the disk of our galaxy. Because it is in the disk rather than in a halo, this dark matter 
must be dissipative which presumably means that it is baryonic. (I am assuming that 
we have necessarily discovered the existence of any form of matter sufficiently abundant 
and sufficiently strongly interacting to be the disk dark matter. Disk matter must 
have sufficiently strong interactions to have concentrated in a disk by dissipating its 
energy while conserving its angular momentum.) Moreover, the recent discovery 

[13] 

of microlensing in the direction of the galactic bulge indicates that our disk has a 
population of low mass (~ 0.1M Q ) compact objects. These may be 'brown dwarfs', i.e. 
stars too low in mass to shine by nuclear burning, which are the most likely hiding place 

for dark baryons J 14 l 

Eq. (1.10) is compatible with the assumption that all the dark matter in galax- 
ies and clusters of galaxies is in baryons. However, we will see in the next section that 
this economical hypothesis runs into difficulties in most scenarios of galaxy formation. 
Eq. (1.10) is not compatible with fig = 1. Hence, if one believes in inflation and in 
standard nucleosynthesis — and both of these are very well motivated — one must con- 
clude that our universe is dominated, at the 90% level, by a form of dark matter which 
is not baryonic. 



2. DARK MATTER CANDIDATES 
2.1. Baryons 

They are known to exist. Moreover, they are known to be a form of dark matter, 
and the nucleosynthesis constraints allow all dark matter associated with galaxies to be 
in baryons. Hence, a conservative hypothesis may be that Q ~ ~ 0.2. (Inflation is 
given up then.) However, that particular scenario has serious difficulties with galaxy 
formation. The point is that the density perturbations in the baryon distribution that 
should produce galaxies cannot grow by gravitational instability till after the epoch 
of recombination at a temperature T rec ~ 4.10^ K. Recombination is when electrons 
combine with ions to form neutral atoms. Before recombination, the baryons are in 
close thermal contact with the photon gas. Because the latter has pressure, the Jean's 
mass is large: 

Mj ~ 1.8- 10 m M Q (n B h 2 r 2 for T > T rec (2.1) 

in this scenario. The Jean's mass sets the critical scale below which density pertur- 
bations do not grow. After recombination, Mj ~ and density perturbations in the 

matter distribution grow on all scales at the rate ~ R where R is the cosmological 
scale factor. The temperature drops as T ~ Galaxies form when ^ ~ 1 on 

r 

the appropriate mass scale, about 10 12 Mq. For this to happen before the present, the 
density perturbations in the baryon distribution at recombination must therefore have 
a minimum amplitude: 



Sp 
P 



Rrec _ Tq_ _ = 3 

rec Ro T rec 4.10 3 K ' V ' ; 
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One should expect accompanying photon temperature fluctuations of the same order of 
magnitude. These would contradict the upper limits < KT 5 ) on the microwave 
background anisotropy. The scenario has additional difficulties due to the diffusion of 
photons, which tends to erase adiabatic fluctuations in the baryon number density, t 1 ^ 

The above model (O ~ ~ 0.2) came first historically and, as we just saw, 
it led to the expectation of large CMBR anisotropies ~ 1CP 3 or 1CP 4 ) which got 
into more and more severe disagreement with the observations. This outcome provided 
a strong impetus for the development of models with cold dark matter (CDM). Indeed, 
the CDM candidates decouple from photons and baryons. As a result, the density 
perturbations in CDM start to grow by gravitational instability as soon as t > t eq where 

t eq is the time of equality of the radiation and matter energy densities (p ra d. ~ R 4 , 

Pmatt. ~ R~ 3 , Prad. = Pmatt. at t eq ). In models where ttcDM is close to one 5 teq comes 
well before t rec - Because there is more time for their growth, the primordial density 
perturbations in these models are smaller than the RHS of Eq. (2.2). 

At any rate, as we saw, some baryonic dark matter is known to exist and there 
may be large amounts of it, up to 0# ~ 0.2. As already mentioned, a likely hiding 

placet 14 ! for these dark baryons is "brown dwarfs", i.e., stars too low in mass to burn 

by nuclear fusion. Paczynskit 1 ^ pointed out that objects of this kind, generically called 
MACHOs for massive compact halo objects, can be searched for by looking for the 
gravitational lensing of background stars by MACHOs that happen to pass close to 

the line of sight. Three collaborations t 13 l have reported compelling candidates for such 

microlensing events. This very exciting development is reviewed by Ansaril 17 ! at this 
meeting. 

2.2. Neutrinos 

Neutrinos decouple in the early universe at a temperature Tp of order a few 
MeV. After that, each neutrino moves freely and hence its momentum decreases with 
the universe's expansion according to: p v ~ Thus, neglecting inhomogeneities, 

the neutrino phase-space density is given by: 

(2n)° e T v (t) + i 

where T v {t) = ^d^tjj and where g v is the number of neutrino spin degrees of freedom. 

In the standard model, each neutrino flavor contributes g v = 2. T v (t) is usually called 
the 'neutrino temperature' although the distribution (2.3) will deviate from a thermal 
one if the neutrino is massive. The photon temperature also decreases according to T 7 ~ 

most of the time. If this were always true, and T v would remain equal. However, 
at a temperature of order 1 MeV, electrons and positrons annihilate, thus reheating the 
photon gas. Because e + e~ annihilation is adiabatic (the process goes back and forth 
very rapidly compared to the Hubble expansion rate), conservation of entropy allows 
one to relate the photon temperature after annihilation to the temperature the photons 
would have had if there had been no annihilation, which is the neutrino temperature. 
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This yields the famous result: 

T v = 3 T 7 (2.4) 

after e + e~ annihilation. Since = 2.73 K today, T v = 1.95 K. This implies in 
particular that the number density of neutrinos and anti-neutrinos today tv+f^o) = 
113/cm^ per neutrino flavor. From this one readily finds that, in extensions of the 
standard model where the neutrinos have small Majorana masses, their contribution to 

the cosmological energy density is! 18 ! 

i 

where the sum is over flavors. Although none of these neutrinos have been observed 
directly or indirectly, we are confident that they are there because the theoretical argu- 
ments for their existence are very simple and conservative. 

So the next assumption we will consider is that neutrinos constitute most of 
the dark matter and hence that they dominate the cosmological energy density. This 
scenario also tends to run into trouble with galaxy formation. As was already mentioned, 
perturbations in the matter density only start to grow after matter-radiation equality. 
In the present scenario, equality occurs when the temperature is of order the neutrino 
mass. Before that the neutrinos are relativistic and their "free-streaming" erases all 
density perturbations in the neutrino fluid on length scales less than the free-streaming 

distance, t 19 l i.e., the distance a typical neutrino travels from the Big Bang till the time 
of equality. The corresponding mass scale is: 

M i , = 4-1O 15 M (— — j , (2.6) 



in 



which is of order the mass in large galactic clusters. The resulting spectrum of primordial 
density perturbations is heavily suppressed on all mass scales less than M v , including 
the mass scale (~ 10 12 Mq) of individual galaxies. If such a spectrum is used as input in 
computer simulations of large scale structure formation, a poor fit to the observations 

results, t 2 ^ The difficulties neutrinos have with large scale structure formation may be 
eased if topological defects, t 21 l such as cosmic strings, are the source of the density 
perturbations because in this case, the density perturbations continue to be created 
after t&q. 

There is however a separate difficulty with neutrinos constituting galactic halosJ 22 ]| 
Liouville's theorem tells us that the phase-space density is constant following the mo- 
tion. Eq. (2.3) therefore implies that the neutrino phase-space density can nowhere 

be larger than \g v - On the other hand, for neutrinos to constitute a galactic halo, 
their velocities must be less than the escape velocity, which for our own galaxy is of 
order 10~^c. The upper limits on the phase-space density and on the velocity imply the 
following upper limit on the physical density: 

-24 



_ 1 10 V ( Vmax \ 3 ( m v \4 
P»,max- 2^3~ {W^J \l9eV ) ' { ' 
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Eq. (1.3), (1.4), (2.5) and (2.7) tell us that neutrinos can only barely be packed tightly 
enough to constitute the Milky Way halo. Dwarf galaxies also have dark matter halos 
but smaller escape velocities. For these galaxies, the neutrino phase-space constraint is 

severely violated J 23 l 



2.3. WIMPs 

WIMPs is an acronym for Weakly Interacting Massive Particles. To be specific, 
consider a massive neutral lepton L. If its mass is less than the temperature 
Td at which it decouples from the thermal bath, then L behaves like a neutrino and 
its cosmological energy density is given by Eq. (2.5) or something very similar to it. 

However, t 24 l if exceeds Td, then for > T > Td, the number density of L 
particles: 

(2^1 Anp dp v /qr^ 



n L (T) = n eq (T) = / Wd V j== (2.8) 



+ 1 

falls off exponentially, as e t . The number density n^(T) tracks its equilibrium value 
n eq (T) as long as the annihilation rate of L particles exceeds the Hubble rate. Thus Td 
is given by: 

WannV}n eq (T D ) ~ H(T D ) (2.9) 
where a an n is the annihilation cross-section of L particles. In Eq. (2.9), the dominant 

- m L 

dependence upon Td is the exponential e t d behaviour of n e g(Tp). As a consequence, 
Td is proportional to up to logarithmic corrections. For cross-sections typical of 
weakly interacting particles, one finds: 

T D ^^m L . (2.10) 
Hence, the cosmological energy density in L particles today: 

PL«o) = m L n L (to) = m L n L (tn) (^V * ™L^\^ (^Y (2-") 

\ Ho J {VannV) A D K 1 D J 

where Nq and Nd are the effective numbers of thermal degrees of freedom today and 
at the decoupling of L particles, conservation of entropy from tp till £q having been 
assumed. Remarkably, the Td dependence on the RHS of Eq. (2.11) cancels out because 

of Eq. (2.10) and because H(T D ) = \J^p(T D ) = sJ^N D ^T^. As a result, the 

contribution of L particles, and more generally that of WIMPs, to the cosmological 
energy density depends almost exclusively upon their annihilation cross-section. One 
finds: 

l2 6-10- 27 cm 3 

n WIMP h 2 ~ — . (2.12) 

{VannV) sec 

For the particular case of a heavy neutral lepton L, one has a a nn ~ G^m^ for < 



mz = 91.2 GeV and a ann ~ for mi > mz- In that case, f^/i 2 ~ -\ in the 



m L 
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range few MeV < < mz with Vt^h? = 1 for ~ 2 GeV, and Q^/i 2 ~ m 2 

when mi > mz with Q^/i 2 = 1 for — 10 TeV. Note that if the WIMP is not its 
own anti-particle, there may be a WIMP-antiWIMP asymmetry, similar to the baryon 
asymmetry. In that case there is an additional contribution to S^wiMp/i 2 , aside from 
the one given by Eq. (2.12). 

The best motivated WIMP candidate is the lightest supersymmetric partner 
(LSP) in supersymmetric extensions of the standard model. Typically, this is a linear 
combination of the photino, the zino and the Higgsino. 

Because WIMPs are non-relativistic from the moment of their decoupling, their 
free-streaming distance is very small and hence their free-streaming does not erase den- 
sity perturbations on any relevant scales. For this reason, WIMPs are called Cold Dark 
Matter (CDM). In contrast, neutrinos, because of their large free-streaming distance, 
are called Hot Dark Matter (HDM). The assumptions of CDM, with &CDM — 1; an d of 
a flat (Zel'dovich-Harrison) spectrum of primordial density perturbations yield a model 

of large scale structure formation^] which has been thoroughly tested and which has 
been, by and large, very successful. However, in light of the COBE measurement of 
the cosmic microwave anisotropy which is also a measurement of the primordial den- 
sity perturbations on the largest scales observable, some modification of the pure CDM 
model may be required. 

If WIMPs constitute the halo of our galaxy, they may be searched for on earth by 
looking for WIMP + nucleus elastic scattering in a laboratory detector. t 26 l The nuclear 
recoil can be put into evidence by low temperature calorimetry, by ionization detection 

or by the detection of ballistic phononsJ 2 ^ WIMPs may also be searched for by looking 
for the decay products (photons, anti-protons. . .) of WIMP annihilation in the halo of 

our galaxyt 2 ^] or by looking for neutrinos produced by the annihilation of WIMPs that 

have been captured by the sunJ 29 ! 



2.4. Axions 

The axion is a hypothetical particle whose existence would insure that the strong 
interactions conserve P and CP in spite of the fact that other interactions violate those 

symmetries. t 3 °] Indeed, the action density of the standard model of elementary particles 
contains in general a term: 

L 9 = ^G^G a ^ (2.13) 

where G^ v is the gluonic field strength, G a ^ v is the dual of Gp U , and g s is the QCD gauge 
coupling. If 9 ^ 0, non-perturbative QCD effects induce violations of P and CP in the 
strong interactions. No such violation has been observed. In particular, the upper limit 
on the neutron electric dipole moment requires 9 < 1CP 9 . But there is no reason in the 
standard model for the parameter 9 to be small. This shortcoming has been called the 
"strong CP problem". 

Peccei and Quinn modified the standard model in such a way that the parameter 
9 in Eq. (2.13) gets replaced by ^j^- where a(x) is a dynamical pseudo-scalar field whose 

J a 

quantum is called the axion ; f a is a quantity with dimension of energy called the axion 
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decay constant. By construction, the vacuum expectation value of a(x) is indifferent 
except for those non-perturbative effects that make QCD depend upon 9. The latter 

produce an effective potential V{9) = V whose minimum is at 9 = 0. Thus 

by postulating an axion, 9 is allowed to relax to zero dynamically and the strong CP 
problem is solved. 

The properties of the axion can be derived using the methods of current algebra. 
The axion mass is related to f a by: 

m a ^ 0.6eV 1Q7GeV . (2.14) 

fa 

All the axion couplings are inversely proportional to f a - Thus, a very light axion is also 
very weakly coupled, t 31 ] A priori, the value of f a: and hence that of m a , is arbitrary. 
However, astrophysical considerations t 32 l and searches for the axion in high-energy and 
nuclear physics experiments^ 33 ! rule out m a > 1CP 3 eV. On the other hand, cosmology 
places a lower limit on m a of order 1CP 6 eV by requiring that axions do not overdose 
the universe. t 34 l 

Indeed, for small masses, axion production in the early universe is dominated by 
a novel mechanism. The point is that the non-perturbative QCD effects that produce 

the effective potential V (-j^T j are strongly suppressed at temperatures high compared 

to Aqcd- At these high temperatures, (a(x)) has arbitrary value. At T ~ 1 GeV, 
the potential V turns on and the axion field starts to oscillate about its CP conserving 
minimum. These oscillations do not dissipate in other forms of energy because, in the 
relevant mass range, the axion is too weakly coupled for that to happen. The oscillations 
of the axion field may be described as a fluid of axions. The typical momentum of these 
axions is the inverse of the correlation length of the axion field at T ~ 1 GeV. Since that 

correlation length is of order the horizon then, we have p n ~ j— ^ — ~ tttt, — ~ 10 - ^ eV. 

to ' ru fiGeV 10~ b sec 

Thus, the axion fluid is very cold compared to the ambient temperature. Its contribution 
to the present cosmological energy density is found to be of order 

^ . 3 (^y(^) f . (2 . 15) 
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Several sources of uncertainty affect the relationship between VL a h? and m a , amongst 
which are the nature of the QCD phase transition and the contribution to O a /i 2 from 
cosmic axion strings. [35,36] Also, if inflation occurs and the post-inflation reheating tem- 
perature is less than f a , then the axion field gets homogenized and there may be an 
accidental suppression of Vt a h? because the axion field happens to lie everywhere close 
to the CP conserving minimum of V. From the point of view of large scale structure 
formation, axions are cold dark matter since they are non-relativistic from the moment 
of their production during the QCD phase transition, as was emphasized above. 

Axion dark matter may be searched for by stimulating the conversion of galactic 

halo axions to photons in a laboratory magnetic field. t 37 l The relevant coupling is 

a a - - , 
L ail = 9l -—E-B (2.16) 
7T fa 
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where is a model-dependent coupling constant of order one. If an electromagnetic 
cavity is permeated by a static approximately homogeneous magnetic field Bq and the 

resonant frequency of the lowest TM (relative to the direction of Bq) mode equals 
the axion mass, some galactic halo axions will convert to quanta of that cavity mode. 
For Bq ~ 10 Tesla and cavity volumes of order lm 3 , the power from these a — > 7 
conversions becomes detectable in a sufficiently short amount of time to allow a search 
over a large axion mass range. Two pilot experiments using this technique have been 

carried outJ 38 ! At present, a second generation experiment^ 9 ! is under construction at 
Lawrence Livermore National Laboratory that will be able to detect dark matter axions 
if their local density is equal to the local halo density given in Eq. (1.4) or higher, if 
g^ > 1 and if m a is in the range 1.3/ueV< m a < 13/ieV. 

3. THE PHASE-SPACE STRUCTURE OF COLD DARK MATTER 
HALOS 

If a signal is found in the cavity detector of galactic halo axions, it will be 
possible to measure the energy spectrum with great precision and resolution because 
all the time that was previously used in searching for the signal can now be used to 
accumulate data. Hence there is good motivation to ask what can be learned about our 
galaxy from analyzing such a signal. 

In many past discussions of dark matter detection on earth, it has been assumed 
that the dark matter particles have an isothermal distribution. Thermalization has been 

argued to be the result of a period of "violent relaxation" t 40 ! following the collapse of the 
protogalaxy. If it is strictly true that the velocity distribution of dark matter particles is 
isothermal, which seems to be a strong assumption, then the only information that can 
be gained from its observation is the corresponding virial velocity and our own velocity 
relative to its standard of rest. If, on the other hand, the thermalization is incomplete, 
a signal in a dark matter detector may yield additional information. 

J.R. Ipser and I discussed^ 41 ] the extent to which the phase-space distribution of 
cold dark matter particles is thermalized in a galactic halo and concluded that there are 
substantial deviations from a thermal distribution in that the highest energy particles 
have discrete values of velocity. There is one velocity peak on earth due to dark matter 
particles falling onto the galaxy for the first time, one peak due to particles falling out 
of the galaxy for the first time, one peak due to particles falling into the galaxy for 
the second time, etc. The peaks due to particles that have fallen in and out of the 
galaxy a large number of times in the past are washed out because of scattering in the 
gravitational wells of stars, globular clusters and large molecular clouds. But the peaks 
due to particles which have fallen in and out of the galaxy only a small number of times 
in the past are not washed out. 

If the fraction of the local dark matter density which is in these velocity peaks 
is sufficiently large, a direct dark matter search, such as the LLNL experiment, may 
be made more sensitive by having it look specifically for velocity peaks. I. Tkachev, 
Y. Wang and I have been studying galactic halo formation with the purpose of obtaining 

estimates of the sizes and locations of the velocity peaks. t 42 l To this end, we have 
generalized the secondary infall model of galactic halo formation to include angular 
momentum of the dark matter particles. This new model is still spherically symmetric 
and it has self-similar solutions. We find that the typical fraction of the local cold dark 
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matter density in any one of the highest energy velocity peaks is several percent. A 
forthcoming paper will give estimates of the highest energy peaks as a function of the 
amount of angular momentum and other model parameters. 
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